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ABSTRACT 



Context. Thermal emission from dust provides a valuable tool for determining important physical properties of the dense structures 

within molecular clouds. 

Aims. We attempt to map the distributions of dust temperature and H 2 column density of IRAS 05399-0121/SMM 1, which is a dense 

double-core system in Orion B9. We also search for substructures within the cores through high-resolution submillimetre imaging. 

Methods. The source was mapped with APEX/SABOCA at 350 fim. We combined these data with our previous LABOCA 870-yum 

data. The spatial resolution of the new SABOCA image, ~ 3 400 AU, is about 2.6 times better than provided by LABOCA, and is 

therefore well-suited to our purposes. We also make use of the Spitzer infrared observations to characterise the star-formation activity 

in the source. 

Results. The filamentary source remains a double-core system on the 3 400 AU scale probed here, and the projected separation 

between IRAS 05399 and SMM 1 is 0.14 pc. The temperature map reveals warm spots towards IRAS 05399 and the southeastern tip 

of the source. Both IRAS 05399 and SMM 1 stand out as peaks in the column density map. A simple analysis suggests that the density 



profile has the form 



,-(2-3j5> 



as determined at the position of SMM 1. The broadband spectral energy distribution of IRAS 05399 



suggests that it is near the Stage 0/1 borderline. A visual inspection of the Spitzer/IRAC images provides hints of a quadrupolar-like 
jet morphology around IRAS 05399, supporting the possibility that it is a binary system. 

Conclusions. The source splitting into two subcores along the long axis can be explained by cylindrical Jeans-type fragmentation, 
but the steepness of the density profile is shallower than what is expected for an isothermal cylinder. The difference between the 
evolutionary stages of IRAS 05399 (protostellar) and SMM 1 (starless) suggests that the former has experienced a phase of rapid 
mass accretion, supported by the very long outflow it drives. The protostellar jet from IRAS 05399 might have influenced the nearby 
core SMM 1. In particular, the temperature map features are likely to be imprints of protostellar or shock heating, while external 
heating could be provided by the nearby high-mass star-forming region NGC 2024. 

Key words. Stars: formation - Stars: protostars - ISM: individual objects: IRAS 05399-0121 - Submillimetre: ISM 



1. Introduction 

Understanding of Galactic star formation is ultimately linked to 
the understanding of the physical properties of dense cores in 
molecular clouds, because they appear to be the direct proge- 
nitors of new stars. It is well established that the formation 
process of low-mass (~ 0.1 — 2 M0) solar-type stars start in 
gravitationally bound starless cores or so-called prestellar cores 
(Ward-Thompson et al. 1994). Because these sources represent 
the physical conditions of dense gas and dust best before the 
influence of star formation, they are especially useful targets 
for studies of the initial stages of the star-formation process. At 
some point, a prestellar core collapses due to its self-gravity, and, 
after a short (~ 10 2 - 10 3 yr) transitional first-hydrostatic core 



stage jLarson 19691 Masuna gaet al. 19 98 , Bate 20111) . a central 
protostar forms in its centre (the so-called second hydrostatic 
core; e.g., Masunaga & Inutsuka 2000). The emission of these 



* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) under programmes 079.F-93 13(A) and 
090.F-9307(A). APEX is a collaboration between the Max-Planck- 
Institut fur Radioastronomie, the European Southern Observatory, and 
the Onsala Space Observatory. 
** Hubble Postdoctoral Fellow. 



young accreting sources is heavily reprocessed by the dense em- 
bedding gas and dust, and consequently they appear as cold, sub- 
mm sources also known as Class objects (Andre et al. 1993. 
2000). 

Based on statistical estimates, the lifetime of the Class 
stage is ~ 1 x 10 5 yr dEvans et al. 20091 Enoch et al. 20091 . 
However, the accretion rate and duration can be highly de- 
pendent on the initial and/or environmental conditions (e.g., 
Vorobyov2010| IQffneretal. 2012bl ). Over time, the protostel- 
lar envelope dissipates through mass accretion onto the forming 
star and ejection of circumstellar material through outflows (e.g., 
Arce & Sargent 2006). The source emission shifts to shorter 
wavelengths, peaking in the infrared, and the object becomes a 
Class I source. The transition from Class to Class I is thought 
to occur when the envelope mass is approximately equal to 
the mass of the central protostar. However, protostellar masses 
cannot be directly measured, and in practice source classifica- 
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tion is very sensitive to the inclination of the protostellar disk 
and outflow (e.g., |Whitney et al. 2003| IRobitaille et"aT2 007 ; 



Offn er et al. 2012bl l. After the envelope is dispersed, further ac 
cretion is minimal, the source is optically revealed, and it enters 
the pre-main sequence (e.g., Wuchterl & Tscha muter 2003) . 

The target source of the present study is a dense core sys- 
tem in the Orion B9 star-forming region that consists of a can- 
didate prestellar core and a Class 0/1 protostar, namely the IRAS 
(Infrared Astronomical Satellite; Neugebauer et al. 1984 1 source 
05399-0121 (hereafter, 105399). This source has been the target 
of several previous studies, and what follows is an overview of 
these studies in chronological order. 

Wouterloot & Walmsley (1986) searched for 22-GHz H z O 
maser emission from 105399 but it was not detected. The source 
was also included in the NH3 survey by Wouterloot et al. (1988), 
where the (1, 1) line was detected at 9.25 km s but only 
an upper limit to the gas kinetic temperature of T^ n < 11.1 
K could be derived. The same authors later included 105399 
as part of their 12 CO/ 13 CO survey of selected IRAS sources 
(I Wouterloot et al. 1989). As can be seen in their Fig. 1, the 
12 CO(l - 0) line shows a blue asymmetric profile with the 
blue peak being 4 K stronger than the red peak. Moreover, the 
13 CO(l - 0) line peaks near the central dip of the double-peaked 
12 CO line. This could be the signature of the infalling gas mo- 
tions (e.g., IZhou 19921 Myers et al. 1996) . However, it is also 
possible that there are just distinct velocity components along 
the line of sight; for example, there is additional CO isotopo- 
logue emission at ~ 2.5 - 3 km s _1 (Woute rloot et al. 1 989; 
their Fig. 1). In addition, the 12 CO line detected by these au- 
thors shows non-Gaussian, high-velocity wing emission of low 
intensity indicative of outflowing gas. 

Lada et al. (1991) performed a CS(2 - 1) survey of the Orion 
B molecular cloud. The CS clump they found to be associated 
with 105399 (No. 30) is commonly called LBS 30 according to 
authors' last names. Within the 5<x emission level, the effective 
radius and mass of LBS 30 were determined to be 0.21 pc and 
85.6 M© (they adopted a distance of 400 pcQ). Further NH 3 stu- 
dies of 105399 were performed by Harju et al. (1991, 1993). The 
NH3(1, 1) map from Harju et al. (1993; their Appendix) shows 
an elongated structure where an additional ammonia condensa- 
tion can be seen southeast of 105399. The average linewidth and 
7\; n they derived for the mapped source are 0.6 + 0.2 km s and 
13.8 + 3.6 K, respectively, and the mass, as estimated from NH3, 
is 14 M© (they adopted a distance of 500 pc). Caselli & Myers 
(1995) mapped the Orion B9 region in the J - 1 -0 line of 13 CO 
and C ls O, and 105399 was found to lie within the C 18 clump 
in the centre of the map (their Fig. 4). 105399 was also included 
in the large 22-GHz H2O maser survey by Codella et al. (2002) 
but, similarly to Wouterloot & Walmsley (1986), no emission 
was detected. To our knowledge, the first millimetre dust con- 
tinuum emission map of LBS 30 was made by Launhardt et al. 
(1996). They used both the SEST 15-m and IRAM 30-m tele- 
scopes to image the 1.3-mm dust emission of LBS 30, and found 
a dust core associated with 105399 and a subcore southeast of 
it, resembling the morphology seen in the NH3 map by Harju 
et al. (1993). Assuming a distance of 400 pc and dust tempera- 
ture of Tdust = 30 K, Launhardt et al. (1996) estimated the total 



mass of the source to be 3.7 M0. LBS 30 was included in the 
multitransition CS study of Lada et al. (1997). Their large veloc- 
ity gradient modelling of the lines yielded a density of 2 x 10 5 
cirT 3 within the region where CS(5 - 4) line was detected, i.e., 
within the effective radius 0.03 pc at d - 400 pc. Bergin et al. 
(1999) mapped the source in the 7=1-0 line of the molecular 
ions N 2 H+, H 13 CO\ and DCO + , and also in C ls O(l - 0) and 
CS(2 - 1) (see their Fig. 1). Interestingly, the line emission of all 
the species peaks towards the southeastern core, rather than to- 
wards 105399. Particularly their N 2 H + map appears very similar 
to that of NH3 ( |Harjuet al. 1993) , which might not be surprising 
given that both species are excellent tracers of the dense inter- 



1 In the present paper, we adopt a distance of 450 pc to the Orion gi- 
ant molecular cloud (Genzel & Stutzki 19891. The actual distance may 
be somewhat smaller as, for example, Menten et al. (2007) determined 
a trigonometric parallax distance of 414 ± 7 pc to the Orion Nebula. 
In this introduction section, however, we quote the distance values (and 
distance-dependent physical properties) used by different authors. 



stellar gas (e.g., Bergin & Langer 1997 1. By utilising the molec- 
ular ion data, Bergin et al. (1999) constrained the fractional ioni- 
sation in 105399 to lie in the range 9.3 xl0~ 8 < x(e) < 1.8 xl0~ 7 , 
with the best-fit value being ~ 1.3 x 10~ 7 . 105399 also appeared 
in the C 18 0/H 13 CO + survey of Aoyama et al. (2001). They des- 
ignated the source as the H I3 CO + clump No. 7, and derived the 
effective clump radius and mass of 0.23 pc and 18 M0 (assum- 
ing d = 400 pc). 

One of the most intriguing features of 105399 is the highly 
collimated HH 92 jet emanating from it, and the associated 
parsec-scale bipolar outflow discovered by Bally et al. (2002). 
This giant outflow, consisting of the Herbig-Haro objects HH 90, 
91, 92, 93, 597, and 598, has a projected length of about 34' (4. 1 
pc at 415 pc). The HH 92 jet was found to exhibit small-scale 
wiggles, and HH 598 at the northwestern end of the outflow was 
found to deviate from the current orientation of the HH 92 jet by 
3° in projection. One explanation for the observed features is that 
105399 is composed of a binary protostar whose orbital period 
modulates the outflow's mass-loss rate and velocity. Connelley 
et al. (2007) observed 2-/im emission features extending 1' to the 
northwest of 105399. These are likely to be the result of shocked 
H 2 emission from the HH 92 jet. 105399 was included in the 22- 
GHz H 2 maser survey of IRAS sources by Sunada et al. (2007) 
but, similarly to the above mentioned studies, no maser emis- 
sion was detected. This is consistent with the fact that among 
low-mass young stellar objects (YSOs), 22-GHz H2O masers, 
which are collisionally pumped in shocked gas, are found to be 
associated with predominantly Class protostars and in excited 
shocks of protostellar jets in the close vicinity of the central star 
(e.g., Furuya2003). 105399 may have already passed through 
the H2O maser phase, but also viewing geometry and variability 
could play a role. 

105399 was designated as the source J054227.9-012003 in 
the SCUBA Legacy Catalogues by Di Francesco et al. (2008). 
The 450- and 850-^m flux densities they reported are 26.06 
and 3.62 Jy, respectively. Ikeda et al. (2009) performed an 
H I3 CO + (l - 0) survey of dense cores in Orion B, and their core 
No. 56 is associated with 105399. The radius and mass they de- 
rived for this H 13 CO + core are 5276 (0.12 pc) and 6.8 M© (at 
d = 470 pc as they assumed). The virial mass they derived, 
5.8 M©, suggests that the H I3 CO + core is near virial equilib- 
rium. More recently, Miettinen et al. (2009; hereafter, Paper I) 
mapped the whole Orion B9 region in the 870-/im dust con- 
tinuum emission with the LABOCA (Large APEX BOlometer 
CAmera) bolometer on APEX. With the peak intensity of 0.81 
Jy beam -1 , 105399 was found to be the second strongest 870-/mi 
source in the region. Similarly to earlier studies, a dense star- 
less core southeast of 105399 was detected. We suggested that 
it is a prestellar core and called it SMM 1 . We also constructed 
the spectral energy distribution (SED) of 105399 by combining 
data from IRAS, Spitzer (24 and 70 p.m), and LABOCA. The 
dust temperature and mass of the cold envelope, and the bolo- 
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metric luminosity were found to be 18.5 ± 0.1 K, 2.8 ± 0.3 
M©, and 21 + 1.2 L©, respectively (d - 450 pc). Based on 
its SED properties, we suggested that 105399 is in a transition 
phase from Class to I, whereas it was previously classified as 
a Class I object (e.g., Bally et al. 2002 and references therein). 
Miettinen et al. (2010; hereafter, Paper II) performed follow-up 
NH3 observations of the Orion B9 cores. The gas kinetic tem- 
perature and one-dimensional non-thermal velocity dispersion 
towards 105399 and SMM 1 were derived to be 13.5 + 1.6 and 
11.9 ± 0.9 K, and 0.24 and 0.27 km s" 1 , respectively. The latter 
values are transonic, i.e., <x NT ~ c s , where c s is the isothermal 
sound speed. We note that the spatial resolution of our NH3 ob- 
servations, 40" or 0.09 pc at 450 pc, corresponds to the typical 
size of dense cores and the sonic scale length of the interstellar 
medium, A s ~ 0. 1 pc, where turbulence transitions from super- 
sonic to subsonic ( |Vazquez-Semadeni et al. 2003 1. Therefore, at 
the sonic scale we expect <x NT ^ c s . 

Employing the new temperature values with the assumption 
that 7\i„ = Tdust, the masses of 105399 and SMM 1 were esti- 
mated to be 6.1 + 1.4 and 8.4 ±1.5 M . The virial-parameter 
analysis of SMM 1 suggested that it is gravitationally bound, 
supporting our earlier speculation that it is in the prestellar phase 
of evolution. Miettinen et al. (2012; hereafter, Paper III) pre- 
sented further molecular-line observations towards 105399 and 
SMM 1. The CO depletion factors towards these sources were 
estimated to be 3.2 + 0.7 and 1.9 + 0.3, respectively. Despite the 
fact that no evidence of a significant CO depletion in SMM 1 
was found, we derived a very high level of N2D + /N2H + deutera- 
tion, namely 0.992 ± 0.267. On the other hand, towards 105399, 
the DCO + /HCO + and N 2 D + /N 2 H + column density ratios were 
found to be 0.020 + 0.004 and 0.207 + 0.046. A lower limit to 
the fractional ionisation of x(e) > 1.5 x 10~ 8 was derived in 
105399, which is over six times lower than the corresponding 
value determined by Bergin et al. (1999). The cosmic -ray ionisa- 
tion rate of H 2 was estimated to be £h 2 ~ 2.6 x 10~ 17 s , which 
is within a factor of two of the standard value 1.3 x 10~ 17 s 
(e.g. JCaselli et al. 19 98). We also detected D 2 CO emission from 
both 105399 and SMM 1. We note that D 2 CO, a molecule ex- 
pected to form on the surface of dust grains, was also detected in 
starless/prestellar cores by Bacmann et al. (2003) and Bergman 
et al. (2011). Its presence in the gas phase requires an efficient 
desorption mechanism, which in starless cores could be due to 
cosmic-ray bombardment and/or the release of the formation en- 
ergy of the newly formed species (see Paper III and references 
therein). However, in SMM 1 the gas-phase D 2 CO could origi- 
nate in shocks driven by the HH 92 jet from 105399. This con- 
forms to the low CO depletion factor derived towards SMM 1 . 
Overall, the results of Paper III suggest very peculiar chemical 
properties in the I05399/SMM 1 -system. 

Recent Herschel data show that I05399/SMM 1 is embed- 
ded in a northeast-southwest oriented filamentary structure (see 
IMiettinen 20121 Fig. 9 therein). Miettinen (2012) found that 
there is a sharp velocity gradient in the region (across the 
short axis of the parent filament), and speculated that it might 
represent a shock front resulting from the feedback from the 
nearby expanding H 11 region NGC 2024. As for the other mem- 
bers of Orion B9, the formation of the I05399/SMM 1 -system 
might have been triggered by this feedback. Finally, we note 
that 105399 is the source No. 3133 in the Spitzer YSO survey 
of Orion molecular clouds (Megeath et al. 2012 ', and No. 310 



shocks (T - 3 548^29 K) near the protostar (at radii > several 
100-1 000 AU). 

To broaden our view of the I05399/SMM 1 -system and its 
physical properties, we mapped it in the 350-yum dust continuum 
emission with the Submillimetre APEX Bolometer Camera 
(SABOCA) on APEX. This allowed us to improve the angu- 
lar resolution by a factor of about 2.6 compared to our previ- 
ous LABOCA map of the source. In this paper, we discuss the 
results of our new high-resolution 350-/im data, and use it in 
conjuction with the LABOCA 870-//m data. The observations, 
data, and data reduction are described in Sect. 2. We present the 
observational results in Sect. 3, whereas analysis and its results 
are presented in Sect. 4. We discuss the results in Sect. 5 and in 
Sect. 6, we summarise the paper and draw our main conclusions. 



2. Observations and data reduction 

2.1. SABOCA 350-p.m imaging 

A field of 7:0 x 6:i (0.012 deg 2 ), centred between 105399 and 
SMM 1 on the coordinates #2000.0 = Q5 h 42 m 29 5 .6, £ 2 qoo,q = 
-01°20'1972, was mapped with SAB OCA (|Siringo et al. 2TJ10) 
on the APEX 12-m telescope (Giisten et al . 20061 1 at Llano de 
Chajnantor in the Atacama desert of the Chilean Andes. This 
site provides very good weather conditions for submm astro- 
nomy (e.g., Tremb lin et al. 20 Oi l. SABOCA is a 37-channel on- 
sky TES (transition edge sensor) bolometer array operating at 
350 yum, with a nominal resolution of 7'.'7 (HPBW). The effec- 
tive field of view of the array is 1!5. The SABOCA passband 
has an equivalent width of about 120 GHz centred on an effec- 
tive frequency of 852 GHz. The observations were carried out on 
22 and 23 August 2012. The atmospheric attenuation at 350 //m 
was monitored using the sky-dip method, and the zenith opacity 
was found to be in the range r z A ™ = 0.84 - 1.17. The amount 
of precipitable water vapour (PWV) was ~ 0.3 - 0.5 mm. The 
telescope focus and pointing were optimised and checked at re- 
gular intervals on the planets Mars and Uranus, the red super- 
giant VY Canis Majoris (VY CMa), FU Orionis star V883 Ori, 
infrared cluster NGC 2071 IR, and the massive YSO/H 2 maser 
source G305.80-0.24 (B13134). The absolute calibration uncer- 
tainty for SABOCA is ~ 30%. Mapping was performed using a 
"fast- scanning" method without chopping the secondary mirror 
dReichertz et al. 20011 >. The mosaic was constructed by combin- 
ing 30 individual fast-scanning (~ 1.'2 - 1'5 s -1 ) maps. The to- 
tal on-source integration time was 4.2 h (during the UTC time 
ranges 08:42-14:02 and 15:26-16:39). 

Data reduction was done with the pipeline iterations of the 
CRUSH-2 (Comprehensive Reductio n Utility for SHARC-2) 
(version 2.12-2) software package^ (Kovacs 2008). During the 
reduction process, smoothing with a Gaussian kernel of the 
size l'/41 was applied. The instrument beam FWHM used in 
CRUSH-2 is 7'.'5, so the angular resolution of the final image 
is 7V63 (0.017 pc or 3 433.5 AU at 450 pc). The gridding was 
performed with a cell size of 1'.'5. The resulting \cr rms noise 
level in the final map is ~ 80 mJy beam 1 . 



2.2. LABOCA 870-^im data 



in the Herschel spectroscopic survey of protostars in Orion 



Our LABOCA ( |Siringo et al. 2009} 870-jum observations were 
already published in Paper I. The data was originally re- 
duced using the BOlometer array data Analysis software 



(Man oj et al. 2013] ). Manoj et al. (2013) observed far-infrared 

CO lines, and suggested that the emission originates in low- 

density gas [n(H 2 ) = 6.3^g2 X 10 3 crrT 3 ] heated by outflow 2 http://www.subituti.caltech.edu/~sharc/crush/index.htrn 
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(BoA; Schulle r20121 l. However, to eliminate the effect of data- 
reduction software, and to be properly comparable with our 
present SABOCA data, we re-reduced the LABOCA data with 
CRUSH-2. Smoothing with a Gaussian kernel of the size 
3'.' 76 was applied, while the instrument beam FWHM used in 
CRUSH-2 is 1975. Therefore, the angular resolution of the fi- 
nal image is 19786 (0.043 pc at 450 pc). The gridding was per- 
formed with a cell size of 470. The resulting l<x rms noise level 
in the final map is 30 mJy beam 1 . For comparison, in Paper I 
the resolution and rms noise were 20713 and 30 mJy beam -1 . 
The resulting peak surface brightnesses of the sources are very 
similar to the values in the BoA-reduced map. 

In Paper III, we noticed that our LABOCA map employed in 
Papers I and II was misaligned, and therefore we had to adjust its 
pointing (see footnote 2 in Paper III). The target positions of our 
previous molecular-line observations were chosen to be the peak 
positions of the LABOCA map before adjusting the pointing, 
and therefore they are slightly offset from the 870-yum maxima. 
The misalignment of the map appears to be the result of our data 
reduction with BoA. In the present work, where the LABOCA 
data is re-reduced with CRUSH-2, there is a good agreement 
between the coordinates of the SABOCA and LABOCA peak 
positions (Sect. 3.2). 

2.3. Archival data from the Spitzer Space Telescope 

In this study, we also employ the infrared data from Spitzer 
dWerner et al. 2 004) to characterise the star-formation activity in 
the I05399/SMM 1 -system. The Spitzer Science Center (SSCfl 
and Infrared Science Archive (IRSAjj recently announced the 
release of a set of Enhanced Imaging Products (SEIP) from the 
Spitzer Heritage Archive (SHAfl These include Super Mosaics 
and a Source List of photometry for compact sources. 

The IRAC (Infrared Array Camera) Super Mosaics were 
downloaded from SHA. The IRAC instrument took simultane- 
ous images at wavelengths of 3.6, 4.5, 5.8, and 8.0 pm with 
an angular resolution between 175 and 179 dFazio et al. 2 004). 
We also retrieved the 24- an d 70-pm MIPS (Multiband 
Imaging Photometer for Spitzer; Rieke et al. 2004i images to- 
wards I05399/SMM 1. The angular resolution at these wave- 
lengths is 6" and 18", respectively. The MIPS data are part 
of the programme "A MIPS Survey of the Orion LI 641 and 
L1630 Molecular Clouds" (PI: G. Fazio; programme name/ID: 
ORION_MIPS/47; AOR-Key: 12647424). 

105399 is designated as SSTSLP J054227.67-012000.4 in 
the Source List. In the 378 diameter aperture, the reported IRAC 
flux densities in mJy are 1.317 ± 0.005 (3.6 pm), 5.915 + 0.001 
(4.5 pm; bandfill value), 7.195 ± 0.006 (5.8 /mi; bandfill), and 
9.796 ± 0.018 (8.0 pm). The 24-^m flux density, from a PSF fit- 
ting method, is 1.153 + 0.524 Jy. The latter is close to the value 
1.3 ± 0.05 Jy determined in Paper I. 

3. Observational results 

3.1. Images 

The SABOCA and LABOCA images are shown in Fig. Q] The 
overall filamentary emission morphology seen in the maps is 
rather similar to each other. The projected length along the long 
axis of the 350-jt/m filament is about 2' or 0.26 pc. The projected 



average width of this filament is about ~ 30" or ~ 0.07 pc. Its 
aspect ratio is therefore roughly A =* 4. The projected separation 
between the 350-/mi peaks of 105399 and SMM 1 is ~ l'l or 
0.14pc. 

The Spitzer images are shown in Figs. [2] and [3] As can be 
seen from the images, 105399 is associated with 24- and 70- 
pm point sources, whereas SMM 1 shows no such emission. 
As described in Sect. 1, 105399 drives the HH 92 protostel- 
lar jet. In the Spitzer/IRAC three-colour image, we have la- 
belled several of the outflow features associated with HH 92 



http : //ssc . spitzer . caltech . edu/ 



http : //irsa . ipac . caltech . edu/ 



http : //sha . ipac . caltech. edu/applications/Spitzer/SHA/ 



(see Bally et al. 2002 Table 1 therein). The very long length of 
the outflow could mean that we are viewing 105399 from an 
edge-on orientation so that the outflow lies in the plane of the 
sky. The 4.5-jum band is particularly useful for studying proto- 
stellar jets as it is sensitive to shock-excited line features (H2 
and CO; see, e.g., jSmith & Rosen 20051 lYbarra & Lada 20091 
De Buizer & Vacca 20 1 01 > . From the base to the northwest tip of 
the HH 92 jet, a chain of several 4.5-//m knots are apparent. The 
brightest H2 knot reported by Bally et al. (2002) is also asso- 
ciated with a strong 4.5-/im feature. These knots might be the 
result of internal shocks within the jet. The knot HH 91 A, show- 
ing extended 4.5-//m emission, is a large shock front, and the 
4.5-//m features of HH 93 show other shocks associated with 
the counterflow (in Fig. [2] we have labelled the closest point to 
105399 and the brightest knot of HH 93). A visual inspection of 
the 3.6- and 4.5-pm images provides hint of a quadmpolar out- 
flow, supporting the suggestion that 105399 is actually a binary 
system. 



3.2. The submm cores 

To extract the 350- and 870-yum core properties, we em- 
ployed the commonly used two-dimensional Clumpfind algo- 
rithm, cl£ind2d, developed by Williams et al. (1994). The al- 
gorithm requires two configuration parameters: i) the intensity 
threshold, i.e., the lowest contour level, which determines the 
minimum emission to be included into the core; and ii) the con- 
tour level spacing, which determines the required "contrast" be- 
tween two cores to be considered as different objects. We set 
both of these parameters to 3<x or 240 mJy beam -1 for SABOCA 
and 90 mJy beam -1 for LABOCA. Blended source structures are 
assigned via a friends-of-friends algorithm to the closest emis- 
sion peak. The stepsize used for contouring the data can affect 
the total number of identified structures: if contour spacings are 
too large real fragments can be missed, while closely spaced 
contours can cause artificial noise features to be identified as 
"real" objects (Pine da et al. 20 09). However, our source is rela- 
tively simple and the Clumpfind results are therefore expected 
to be robust. Moreover, Pineda et al. (2009) concluded that that 
the 2D version of Clumpfind used here is more reliable than the 
three-dimensional version. 

The J2000.0 coordinates of the peak emission, angular and 
linear core effective radii (R e g = ^jA pm] /n, where A pro j is the 
projected area within the 3cr contour), peak surface brightnesses, 
and integrated flux densities (within 3cr) are listed in Cols. (3)- 
(8) of Table Q] The effective radii listed in Table Q] are not 
corrected for beam size. The quoted flux density uncertainties 
were derived from the rms noise and the absolute calibration er- 
ror (30% for SABOCA and 10% for LABOCA), which were 
combined quadratically. We would like to point out that the 
LABOCA emission peak of 105399 lies 878 southeast of the pre- 
viously determined position (in the BoA-reduced map), whereas 
for SMM 1 the revised maximum lies 6" to the east. 
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Fig. 1. SABOCA 350-yum and LABOCA 870-//m images of the I05399/SMM 1 -system. The images are shown with linear scaling, 
and the colour bars indicate the surface-brightness scale in Jy beam 1 . The overlaid SABOCA contours start at three times the 
noise level (3<x = 240 mJy beam 1 ) and increase at this interval to 1 680 mJy beam 1 . The LABOCA contour levels range from 
90 mJy beam" 1 (3<x) to 720 mJy beam -1 , increasing in steps of 3cr. The green plus signs show the target positions of our previous 
molecular-line observations. A scale bar indicating the 0. 1 pc projected length is shown in the left panel, with the assumption of a 
450 pc line-of-sight distance. The effective beam sizes (7'.'63 for SABOCA, 19'.'86 for LABOCA) are shown in the lower left corner 
of the panels. 



Table 1. SABOCA 350- and LABOCA 870-yum characteristics of 105399 and SMM 1. 



Source 


A 


a 2000.0 


^2000.0 


R 


eff 


.peak 


Sv 




[jim] 


[h:m:s] 


[":':"] 


["] 


[pc] 


[Jy bearrT 1 ] 


[Jy] 


IRAS 05399-0121 


350 


05 42 27.7 


-01 20 00.4 


19.3 


0.042 


6.24 ± 1.87 


20.66 ± 6.20 




870 


05 42 27.7 


-01 19 57.0 


35.3 


0.077 


0.82 ± 0.09 


2.09 ± 0.21 


SMM1 


350 


05 42 30.9 


-01 20 45.4 


18.7 


0.041 


1.03 ± 0.32 


7.96 ± 2.39 




870 


05 42 30.9 


-01 20 45.0 


33.1 


0.072 


0.40 ±0.05 


1.49 ±0.15 



4. Analysis and results 

4.1. Dust temperature and H 2 column density maps 

The present 350-//m data together with our previous 870-yum 
data allow us to determine the dust colour temperature ( Tdust) and 
H 2 column density [iV(H 2 )] distributions in the I05399/SMM 1 

-system. 



We smoothed the SABOCA map to the resolution of our 
LABOCA data (7/63 -> 19'.'86), and regridded the maps 
to the same pixel size (4"). Convolving the SABOCA map 
halved the \cr rms noise level to 40 mJy beam -1 . The inten- 
sity ratio map was then converted into a Tjust map with the 
dust emissivity index fixed to j3 - 1.8. The choice of this 
value of yS, which lies between the usual fiducial values of 
j3 — 1 - 2, is explained below. Fitting fi would require more 
than two photometric bands, particularly on the Rayleigh- Jeans 
side of the SED (e.g., |Doty & Leung 1994| |Shetty et al. 2009l >. 
Unfortunately, the SCUBA Legacy 450-jum map does not cover 
the I05399/SMM 1 -system completely (about half of SMM 1 
lies outside the map's boundaries). The Orion B9 region was 
observed as part of the "Herschel Gould Belt Survey (HGBS)" 



(I Andre et a l. 20100 Besides the value of /? and the paucity of 
wavelength bands used, another important factor affecting the 
derived dust colour temperatures is that our 870-yum continuum 
data could be significantly contaminated by the 12 CO(3 - 2) line 
emission at 345 GHz (Dra beket al. 20121 1. Although CO does 
not appear to be significantly depleted in the studied source 



(Bergin et al. 1999 Paper III), the CO contamination is not ex- 
pected to be more than a few tens of percent. Mapping obser- 



6 Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA (Pilbratt et al. 2010). 
The HGBS is a Herschel Key Programme jointly carried out by 
SPIRE Specialist Astronomy Group 3 (SAG 3), scientists of sev- 
eral institutes in the PACS Consortium (CEA Saclay, INAF-IFSI 
Rome and INAF-Arcetri, KU Leuven, MPIA Heidelberg), and sci- 
entists of the Herschel Science Center (HSC). For more details, see 
http : //gouldbelt-herschel . cea . fr At the moment, the unpub- 
lished Level-2 map products are available in the Herschel Science 
Archive (HSAjJ but these data do not have an absolute intensity level 
(M. Hennemann, private communication). Unlike Herschel, the ground- 
based bolometers SABOCA and LABOCA heavily filter the emission 
on large scales in order to remove the atmospheric emission. Because 
of the current quality of the Herschel data on I05399/SMM 1, we do not 
employ them in the present study. 
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vations of CO(3 - 2) would be required to quantify this effect. 
The constructed map of line-of-sight-averaged dust temperature 
is shown in the left panel of Fig. |4] 

As expected, 105399 appears as a warm spot with r^st val- 
ues up to ~ 50 K in the temperature map due to embedded heat- 
ing source. A higher /3 value would lower the temperature for the 
same flux density ratio. Within the aperture with size equal to the 
effective 350-/mi angular size (R e s - 19'.'3), and centred on the 
SABOCA peak, the average r dus t towards 105399 is 22.3 ± 10.2 
K, where the + -error represents the standard deviation (hereafter 
referred to as std). This is comparable to the value ~ 18.5 K ob- 
tained from the SED fit in Paper I. On the other hand, the mean 
and std of T^st is 15.3 + 7.0 K in a 40"-diameter aperture centred 
on the position of our previous NH3 measurements with 40" res- 
olution (Paper II). This is close to the measured gas temperature 
of 13.5 + 1.6 K. As can be seen in Fig.|4l the temperature is colder 
in other parts of the system. Measuring r^st within R e ff = 18" 7 
and centred on the SABOCA peak of SMM 1, the mean+std is 
12.5 ± 1.2 K. Towards our line observation position, the corre- 
sponding value is 11.4 + 3.0 K within a 40" aperture. This is 
strikingly similar to the 7\i n value of 1 1 .9 ± 0.9 K derived from 
ammonia. The T^ust rises to ~ 15-16Kat the southeastern edge 
of the source. 

The ./V(H 2 ) map was constructed as follows. We first com- 
puted the optical thickness map at 870 //m. The t^o,™ map 
was then divided by the 870-jum dust extinction cross-section 
per H 2 molecule [cf. Eq. (9) in Miettinen & Harju (2010)]. 
We adopted the solar abundances for hydrogen, helium, and 
heavier elements, i.e., their mass fractions were taken to be 
0.71, 0.27, and 0.02 of the total mass of all elements, re- 
spectively. This corresponds to a helium to hydrogen abun- 
dance ratio of about 0.10 and mean molecular weight per H2 
molecule of p^ - 2.82 (Kauffm ann et al. 2 008; Appendix A.l 
therein). The 870-/im dust mass absorption (or emission) co- 
efficient, i.e., dust opacity per unit dust mass was taken to be 
1.38 cm 2 g -1 . This value was interpolated from the widely used 
Ossenkopf & Henning (1994, hereafter OH94) model describ- 
ing graphite-silicate dust grains that have coagulated and ac- 
creted thit\2 ice mantles over a period of 10 5 yr at a gas den- 
sity of « H = «(H) + 2n(H 2 ) =* 2n(H 2 ) = 10 5 cm 4 . For the 
average dust-to-hydrogen mass ratio, M^ ufil /Mu, we adopted the 
canonical value 1/100 (e.g., Drai ne 20111 Table 23.1 therein), al- 
though the true value can be different (e.g., Vuong et al. 2003 
IDraine eFa l. 2007). The total dust-to-gas mass ratio used is 
therefore M dust /M gas = M d ust/(1-41M H ) = 1/141, where the fac- 
tor 1.41 (^ 1/0.71) is the ratio of total mass (H+He+metals) to 
hydrogen mass. In the adopted OH94 dust model, /? =* 1 .8, as de- 
termined from the slope between 350 and 1 300 /urn (ka °c A~@). 
Therefore, this value was adopted in the derivation of the rd ust 
map to be consistent with the dust properties used. The result- 
ing column density map is shown in the right panel of Fig. |4] 
Both 105399 and SMM 1 stand out as dense portions of the fil- 
amentary structure. The 2V(H 2 ) values towards these two cores 
are ~ 5 - 6 x 10 22 and ~ 3 - 4 x 10 22 crrT 2 , respectively. 



4.2. Distribution of the dust emissivity index 

We can also estimate the value of the dust emissivity index 
by employing a constant value for the dust temperature (e.g., 



Spitzer/IRAC 3.6 |im 



WJy/sr 



8 In Papers I — III, we assumed that the grains have thick ice mantles, 
and used the 870-yum opacity of K87(Vm — 1.7 cm 2 g _1 from OH94. The 
assumption of thin ice mantles is supported by the fact that no signifi- 
cant CO depletion was found in Paper III. 
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Fig. 3. A zoomed in view of Fig. [2] to show the Spitzer/IRAC 
3.6-yum emission towards 105399. The image is displayed with 
logarithmic scaling to highlight the jet features. The dashed lines 
illustrate the possible quadrupolar outflow structure. Note that 
SMM 1 lies outside the image. The black plus sign shows the 
peak of the 24-jum source. 



ISchnee & Goo dman 2005]). As mentioned in Sect. 1., Harju 
et al. (1993) found that the average gas temperature in the 
I05399/SMM 1 -system is 13.8 + 3.6 K, without significant vari- 
ation between different parts within the errors (J. Harju, private 
communication). Indeed, the 7^ values determined in Paper II 
towards 105399 and SMM 1 are both close to this average value. 
Therefore, it seems rather reasonable to assume that the aver- 
age dust temperature of the system is about 13.8 K. The derived 
emissivity spectral index map is shown in Fig. [5] 

Measuring again the values towards the 350-^m peaks with 
aperture sizes corresponding to the core effective sizes, we ob- 
tain the average j3 values of 2.2 + 0.4 and 1.6 + 0.2 for 105399 
and SMM 1, respectively. These values are physically reason- 
able, and, within the errors, comparable to value /? = 1.8 adopted 
above. 



4.3. Core masses, and H 2 column and volume-averaged 
number densities 

The core masses, M, over an effective area of radius /? e ff were es- 
timated from the integrated 350- and 870-yum flux densities using 
the standard optically thin dust emission formulation [see, e.g., 
Eq. (2) in Paper I]. As the dust temperature we used the average 
values 22.3 + 10.2 K and 12.5 ± 1.2 K for 105399 and SMM 1, 
respectively. Following the OH94 dust model described above 
(thin ice mantles), the dust opacity per unit dust mass at 350 and 
870 /vm was taken to be K/^o^m = 7.84 cm 2 g~' and/s^o^m = 1.38 



J/jm 



cm~ g , The total dust-to-gas mass ratio was taken to be 1/141. 
The uncertainty in mass was propagated from the uncertainty in 
flux density, and for SMM 1 also from that of the temperature 
(the large std of T^st for 105399 was not taken into account). 
The uncertainty in dust opacity, which is likely being a factor 
> 2, is the major source of error in the mass estimate. 
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Fig. 5. Distribution of the dust emissivity spectral index with the 
assumption that the dust temperature is constant at 13.8 K (= 
(7\i„); Harju et al. 1993). The diamond symbol indicates the 24- 
yum emission peak of 105399. 



lArzoumanian et al. 201 11 cf. Sect. 4.5). In the case of a sim- 
ple spherical source with a power-law volume (number) den- 
sity profile of the form n(H2) <* r~ p , the column density is ex- 
pected to follow the relation N(H2) °c r~ m with m - p - 1 (e.g., 
lYun & Clemens 199 It . In our case, the density profile would 
then be ra(H.2) °c r - ' 2 ' 3 ^). The deconvolved size of SMM 1, 
when determined from the LABOCA image [whose resolution 
equals that of the N(Ha) map], is very close to the beam size 
19'.'9. Interestingly, when this fact is compared with the density 
power-law distribution models from Young et al. (2003; their 
Fig. 27), a value of p =* 1.9 would be expected, not far from the 
value deduced above (cf. Sect. 5.3 in Paper I). 



5.0x10^ 



3.8x10" - 



E 2.8x10 - 



2.1 x10^ 




Besides the values extracted from the A^PL) map above, 
we also computed the peak beam-averaged H2 column densities 
from the peak surface brightnesses in a standard way [see, e.g., 
Eq. (3) in Paper I]. The parameters needed in the calculation 
(Tdust, Ph 2 , k a, A^dust/A^gas) were the same as described above. 
The error in N(H2) is based on the uncertainty in the peak sur- 
face brightness, and for SMM 1 also on the 7d us t uncertainty. 
Similarly to mass calculation, the k,\ uncertainty is likely to be 
the main source of error here. 

The volume of a sphere-averaged H2 number densities over 
R e fi, (n(H2)), were calculated using Eq. (1) of Paper III, and 
the corresponding errors were propagated from those of M. The 
masses and densities derived above are listed in Cols. (4)-(6) of 
Table [2] The core masses derived from 350-/mi data are about 
half the values derived from the lower resolution 870-yum data. 
On the other hand, the column and number densities derived 
from 350-yum data are typically about three times higher com- 
pared to those derived from LABOCA data. 

We also constructed a column density distribution us- 
ing azimuthally averaged values (see Fig. |6j. As shown in 
Fig. |U the 4"- or one pixel-wide annuli were concentric cir- 
cles centred on the A^PL) peak of SMM 1 at (a, 5)20000 = 
5 h 42 m 30 s 37, -1°20'32'.'95. The radial N(H 2 ) distribution is pre- 
sented in Fig. [6] In order to characterise the structure, the 
data points were fittecQ with the function [cf. Eq. (15) in 
Nielbock e t al. 201211 



JV(H 2 ) 



JV(H 2 )o 



1 



(t) 



(i) 



where A^PDo is the central peak column density, r is radial dis- 
tance, ro is the radius of the flat inner region, and m is the power- 
law exponent at large radii (r » ro). The obtained values are 

N(H 2 )o = (4.1!};J) x 10 22 cirT 2 , r = 28'.'3+^ 6 * 0.06!°$ 

pc, and m - 1.3^H. We note that ro is only roughly compa- 
rable to the thermal Jeans length at the centre of the annuli, 
i.e., A] - 0.01 ± 0.004 pc at T = 11.9 ± 0.9 K (see, e.g., 



9 We used the non-linear least squares IDL fitting routine MPFIT 
( Markwardt 2009 1 available at http : //purl . com/net/mpf it 
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Fig. 6. Radial profile of the column density of SMM 1 displayed 
with logarithmic scale for both axes. The data points represent 
azimuthally averaged values inside 4" wide annuli (Fig. [4]). The 
vertical error bars reflect the standard deviation of the azimuthal 
averaging. The mean radial A^Fb) profile can be fitted by a 
power-law with A^PL) °c r~(' 3 -o«) between the flat inner core 
and the "edge", as indicated by the thick solid line. 



4.4. Modelling the spectral energy distribution of IRAS 
05399-0121 

To characterise the physical properties of 105399, we constructed 
its SED by combining the Spitzer photometric data with the 
submm data from SABOCA, LABOCA, and SCUBA. We also 
employed the IRAS flux densities (see Table 1 in Paper I), for 
which a typical uncertainty of 10% was adopted. Di Francesco 
et al. (2008) reported a rather large 850-yum effective radius 
for 105399, namely 49V4. The reported SCUBA flux densities 
are therefore likely to include contributions from SMM 1 also. 
For this reason, we used apertures of the size corresponding to 
the effective 350- and 870-/im core sizes to extract the 450- 
and 850-/im flux densities, respectively. The obtained values are 
1 1 .06 + 5.53 and 2.59 ± 0.52 Jy, where the quoted uncertainties 
represent the 50% and 20% absolute flux uncertainties, respec- 
tively. 

To build the SED, we used the SED fitting tool developed 
by Robitaille et al. (2006, 20070 The tool employes a grid 
of two-dimensional axisymmetric radiative-transfer models of 
YSOs ( |Whitney et al. 2003) > to fit the input flux densities. The 



10 The SED fitter is publicly available 

http : //caravan . astro . wise . edu/protostars/ 



at 
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Table 2. Dust colour temperatures, masses, beam-averaged peak H2 column densities, and volume-averaged H2 number densities 
of the cores. 



Source 


(Tdust) 


A 


M 


JV(H 2 ) 


<n(H 2 )> 




[K] 


ljum] 


[M Q ] 


[10 22 cm- 2 ] 


[10 4 cm" 3 ] 


IRAS 05399-0121 


22.3 ± 10.2 


350 


2.1 ±0.6 


8.9 ±2.7 


12.9 ±3.7 






870 


3.8 + 0.4 


3.1 ±0.3 


3.8 ± 0.4 


SMM 1 


12.5 ± 1.2 


350 


3.8 ±1.7 


7.0 ±3.2 


25.2 ±11.3 






870 


6.7 ±1.3 


3.8 ±0.8 


8.2 ±1.6 



grid encompasses a wide range of stellar masses (0.1-50 Mq; 
mainly covering low-mass YSOs) and YSO evolutionary stages 
(stellar ages of 10 3 - 10 7 yr), and makes use of 20000 sets of 
14 physical parameters and 10 different viewing angles (inclina- 
tions ranging from edge-on to pole-on or from 90° to 0°), result- 
ing in a total of 200000 model SEDs. The models assume that 
the central embedded source is forming as a single star, it is as- 
sociated with a flared disk, and that the system is surrounded by 
a dusty envelope with cavities carved by bipolar outflows. Here, 
a dust-to-gas ratio of 1/100 is assumed. The models take into ac- 
count the fact that silicate dust can produce the 10- and 18-/im 
features (stretching of the Si-0 bonds and bending of O-Si-O, 
respectively), but the models do not include the icy grain man- 
tles or PAH emission features at mid-infrared. The model fitter 
requires the visual extinction and distance ranges within which 
to fit the observational data. We limited the former to be within 
the range Ay = ... 8.5 mag, which is based on the median Ay 
value towards the NGC 2024 cluster (IHaisch et al. 2001b . and 
the distance was restricted to 450 pc . 

The SED is presented in Fig. [7] The best fit is shown by a 
black curve, whereas the grey curve shows the next-best fit sat- 
isfying the criterion^ 2 -^ e8t < 5Af d at a , where 7Y data = 13 is 
the number of flux data points. The fit to the stellar photosphere 
SED is shown by a dashed line. Due to possible contamination 
by shock emission, the Spitzer/IRAC 4.5-yum flux density was 
finally neglected from the fit (it was treated as an upper limit 
with a zero weight). Inclusion of this data point raised the lumi- 
nosity to an unrealistically high level (~ 230 L0). Some other 
IRAC bands could also be influenced by PAH emission or scat- 
tering: the 3.6, 5.8, and 8.0 //m bands contain PAH features at 
3.3, 6.2, and 7.7 and 8.6 fim, respectively (Draine 2003), and 
the 3.6-fj.m emission can originate from scattered interstellar ra- 
diation dSteinackeretal.2010l i. The IRAS 60- and 100-/mi flux 
densities appear to be slightly inconsistent with the Spitzer 70- 
yum flux density measured by PSF fitting in Paper I. This is likely 
to result from the fact that the IRAS beam of a few arcminutes 
picks up emission from the source surroundings. 

The corresponding SED parameters are listed in Table [3] In 
this table, we give the inclination angle (/), visual extinction be- 
tween the observer and the outer edge of the circumstellar ma- 
terial (Ay), stellar age, mass, and temperature (t*, M*, T*), the 
rate of mass accretion from the envelope onto the disk (M env ), 
disk mass (M^k), visual extinction of the circumstellar material 
(Ay rc ), total luminosity (L tot ; contributions from the central star 
and accretion, assuming that all the accretion energy is radiated 
away), envelope mass (M env ), and the implied evolutionary stage 
of the source. We give the best-fit model parameters with errors 
calculated as difference between the minimum and maximum 
values that can be derived from the three next-best model fits. 
We note that all the masses and the mass accretion rate assume 
a dust-to-gas ratio of 1/100. The source inclination, 31°8 , 
is roughly consistent with the observed outflow orientation. The 
best-fit model stellar mass and total luminosity are found to be 



M+ ~ 0.5 Mq and L tol ~ 15.5 Lq. The latter is corrected for 
foreground extinction and does not depend on the viewing an- 
gle. It is also comparable to the value 21 + 1.2 Lq derived by 
fitting a two-temperature component modified blackbody to the 



source SED (Paper I). The envelope mass, M el 



1.5 



o.') 



Mq, 



determined from the long wavelength emission, is very similar to 
that derived from the 350-yt/m flux density alone, especially when 
the factor of 1.41 difference in the dust-to-gas ratio is taken into 
account. It is also comparable to the value 2.8 + 0.3 Mq derived 
for the cold SED component in Paper I (where the dust-to-gas 
ratio 1/100 was adopted and grains were assumed to have thick 
ice mantles). 

In the case the range of values is very large, the parameter in 
question is not well constrained (such as the disk mass). 105399 
can be classified following the YSO classification scheme of 
Robitaille et al. (2006). The evolutionary stage given in the last 
row of Table [3] is based on the values of M env and Mdi s k rela- 
tive to M*. 105399 has M en v/M* * 1.4 x 10~ 4 > lO" 6 yr 1 , 
and is therefore a Stage I candidate (dominated by a large ac- 
creting envelope). This stage is taken to encompass both the 
Stage (i.e., 30-K modified black-body SED£J and I sources 
(Robitaill e et al. 20061) . The earliest evolutionary stage is cha- 
racterised by an envelope whose mass is higher than those of 
the central star and the disk (see McKee & Offner 2010). For 
105399 we obtain the M em /M+ ratio of only ~ 3, supporting the 
idea that it is at the Class (or Stage) 0/1 borderline. Offner et al. 
(2012b) simulated the formation of low-mass stars including the 
effects of outflows and radiative heating. They found that while 
the SED models could correctly identify the evolutionary stage 
of the synthetic embedded protostars, the disk and stellar param- 
eters inferred from the SEDs could be very different from the 
actual simulated values. In particular, the fact that 105399 may 
be a binary system makes many of the derived source properties 
highly uncertain. 

4.5. Fragmentation of the source 

To investigate the dynamical state of the detected 350-//m fila- 
mentary structure as a whole, we estimated its total mass as 



M = jt/H 2 m H J] [N(H 2 )A pix ]. 



(2) 



where m^ is the mass of a hydrogen atom, A plx is the surface 
area in one pixel (16 □"), and the sum is over all pixels within 
the 3cr contour of 350-/im emission. This integration of the mass 
surface density gives M ^ 8.8 Mq. As the projected length of 
the source along the long axis is 0.26 pc, the corresponding mass 
per unit length, or line mass, is Mn ne =* 34 Mq pc . 



11 The Robitaille et al. (2006) SED models do not take into account 
the cold envelope material (< 30 K), which might have some effect on 
the longer wavelength fits (see also Offner et al. 2012b). Therefore, the 
mass of the coldest envelope parts are neglected. 
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Fig. 7. Spectral energy distribution (SED) of 105399. The filled 
black circles associated with vertical error bars represent the 
Spitzer, IRAS, SABOCA, SCUBA, and LABOCA flux densities. 
The downwards pointing triangle indicates the Spitzer/IRAC 
4.5-jum flux density considered as an upper limit (not taken into 
account in the fit). The solid black line indicates the model that 
gives the best fit to the input flux densities, whereas the grey 
line illustrates the subsequent good fit model with the criterion 
X 1 -;fb e st < 5» where x 2 is taken per data point (A^ata = 13). The 
black dashed curve shows the SED of the stellar photosphere 
corresponding to the central source of the best-fitting model (as 
it would appear in the absence of circumstellar dust, but includ- 
ing the effects of interstellar dust). Note the presence of silicate 
absorption features at 10 and 18 yum. 

Table 3. Results of SED modelling of IRAS 05399-0121. 



Parameter 


Value" 


Min.^/W data 


103.6+ 36 ' 3 


Inclination, i [°] 


31.8+ 1 " 


Visual extinction, A\ [mag] 


2 2 +2 - 2 

^ -1.4 


Stellar age, T* [10 3 yr] 


6 3 +84 


Stellar mass, M* [Mq] 


0.46+' 09 


Stellar temperature, T+ [10 3 K] 


3.65+ 053 


Envelope accretion rate, M e „v [10~ 5 Mq yr _1 ] 


6 47+9.13 
"■ H '-4.35 


Disk mass, M^ [Mq] 


0352+ 01278 
U - UJJ -0.0345 


Circumstellar extinction, A?) lc [mag] 


430+g 6 


Total luminosity, L to i [L©] 


15 5 +23 - 6 

J - -8.62 


Envelope mass, M em [Mq] 


1 51 +2 - 07 

'• J1 -0.9I 


Stage 


I 



Notes. <a) For each model parameter, the best-fit value is given with the 
± errors computed from the minimum and maximum values obtained 
from the three next-best model fits. For some parameters, the best-fit 
value is the lowest one. 



As having a filamentary shape, the I05399/SMM 1 -system 
can be considered a cylindrically shaped object. It is therefore 
interesting to examine if the source's substructure can be un- 
derstood in the context of cylindrical fragmentation. For an infi- 
nite, unmagnetised, isothermal cylinder, the instability (collapse 
to line singularity) is reached if its Mr lne exceeds the critical equi- 
librium value of (e.g., Ostrik e~r~1964l Inutsuka & Miyama 1992) 



^"e 



2c 



(3) 



where G is the gravitational constant. Using the average gas 
temperature of <7/ kin > = 13.8 ± 3.6 K ( |Harju et al.1993"] ) and 
the mean molecular weight per free particle of p p = 2.37 
(IKauf fmann et al. 2008 (3 to calculate the sound speed, we de- 
rive the value MJ* a 22 + 6 M© pc 1 . The I05399/SMM 1 
-system therefore appears to be a thermally supercritical fila- 
ment by a factor of 1.5 + 0.4. If the external pressure is not 
negligible, as may well be the case for the dynamic Orion 
B9 region, then MSj is smaller than the above value (e.g., 
Fiege & Pudritz 2000). On the other hand, if the non-thermal 
motions are taken into account, i.e., c s is replaced by the effec- 
tive sound speed c e ff = (c 2 + cr^ T )^ 2 , where <x NT ~ c s (Paper 
II), the value of My"l becomes about two times larger than in 
Eq. © ( |Fiege & Pudritz 2000). We would like to stress that this 
includes the assumption that all of the non-thermal motions are 
providing additional support, similarly to that provided by ther- 
mal pressure (IKirket al. 20 13k The fact that our source already 
contains a protostellar object, i.e., is collapsing, suggests that it 
should be (slightly) supercritical. We note that the source mass 
determined from the A^FF.) map suffers from uncertainty in the 
dust properties (such as the dust-to-gas ratio and dust opacity). 

As mentioned in Sect. 3.1, the projected separation between 
the subcores 105399 and SMM 1 is 0.14 pc at d = 450 pc (the 
spacing ranges from 0.13 to 0.16 pc for the distance range 400- 
500 pc). In the case the cylinder has My me - Mj?*, as is roughly 
the case for our source, the Jeans length along the long axis is 
(ILarson 1985llHartmann 20021 ) 



4 = 



3.94c;: 



(4) 



where So = (Jn 2 m nN(fl2) is the central surface density. We esti- 
mated the central column density to be the average value along 
the slice shown in Fig. [4] (position angle 135°6 east of north), 
i.e., ~ 2.6 + 1.0 x 10 22 cirT 2 , where the uncertainty is given 
by one standard deviation. The corresponding value of A c , is 
0.075 ± 0.035 pc. This is about 1.9 + 0.9 times smaller than 
the observed value but within the errors they are in reasonable 
agreement. The system may therefore lie close to the plane of 
the sky. For comparison, the traditional thermal Jeans length is 



A,= 



\—L- 
'G<p>' 



(5) 



where (p) = pu 2 mu(n(H2)) is the mean mass density (e.g, 
Kauffman n et al. 20101) . If we use the average density of ~ 6 ± 
1 x 10 4 crrT 3 derived from the 870-/mi values listed in Col. (6) 
of Table 12 we obtain Aj =* 0.09 ± 0.01 pc. This is also roughly 
comparable to the observed core separation. 

If the finite length and diameter of the cylindrical cloud 
are L and D, respectively, the number of fragments form- 
ing due to gravitational instability is expected to be (e.g., 
Basti en et al. TWTllWiseman & Ho 19981 



AW = T— 



2 (L/D) 



0.5A. 



(6) 



Ac* 3.94 

where A cr [i - 1 .97D is the wavelength of the most unstable per- 
turbation, and A the cylinder's aspect ratio. As mentioned earlier, 



12 The classical value of // p = 2.33 applies for a gas consisting of H 2 
molecules with 10% He (He/H= 0.1) and a negligible amount of metals. 



Miettinen & Offner: SABOCA and LABOCA imaging of IRAS 05399-0121 



A is about 4, so one would expect to have two subfragments, in 
agreement with the two subcores 105399 and SMM 1. Note that 
the average width of the SABOCA filament, D ~ 0.07 pc, yields 
^crit — 0.14 pc, which is in excellent agreement with the ob- 
served separation distance. 

To conclude, the studied system has likely fragmented 
into subcores as a result of cylindrical gravitational instability. 
Considering a filament of radius R, the fragmentation timescale 
is then expected to be comparable to the radial crossing time, 
Across = R/c, where <x refers to the total velocity dispersion. 
If we use the average FWHM linewidth of Av = V81n2cr = 
0.6 km s from Harju et al. (1993), which is very similar to 
the NH3(1, 1) linewidths found in Paper II, the T closs for the 
SABOCA filament is ~ 1.34 x 10 5 yr. If 105399 has accreted 
mass at a constant rate since the formation of the central proto- 
star, the fragmentation timescale is 21*J2 ] ti mes longer than the 
protostar's age inferred from SED modelling (Table |3}- 

5. Discussion 

5.1. Dust properties 

As explained in Sect. 4.1, the temperatures measured for the gas 
component in Paper II are similar to the dust temperatures mea- 
sured here from the 350-to-870-/mi flux density ratios. The sim- 
ilarity between 7\i n and Tdust suggests that the gas and dust are 
well-coupled. This is not necessarily surprising, because ther- 
mal coupling between the gas and dust by collisions is expected 
at densities n > 10 4 cirT 3 (Goldsmith & Langer 1978] >. 

The constructed T& usl map shown in Fig. [4] reveals a warm 
spot towards 105399. Given the presence of an embedded proto- 
star, the increased 350/870-yum flux density ratio is likely domi- 
nated by temperature, although it could be partly due to variation 
in f3 also. However, the 24-fj.m peak position and the warmest 
spot are not exactly coincident, but there is a slight offset of 
about 9" between the two. This difference could be caused by 
the outflow feedback from 105399, namely heating by outflow 
shocks. Also, the Tdust values appear to rise when going south- 
east of the 350-//m peak of SMM 1 . Perhaps this is another re- 
flection of the outflow interaction related to the HH 92 jet. This 
conforms to the low CO depletion factors and to the detection 
of deuterated formaldehyde (D2CO; a species forming on icy 
grain surfaces) reported in Paper III as described in Sect. 1 . It 
is worth noting, however, that the warmer southern parts of the 
I05399/SMM 1 -system face the direction towards the Flame 
Nebula NGC 2024, an active massive star-forming region (see 
Fig. 9 in Miettinen 2012). Even though this could enhance the 
strength of the local interstellar radiation field and thus provide 
external dust heating, one would then expect to see a more uni- 
form Tdust gradient across the short axis and along the whole 
source. Therefore, protostellar or shock heating seems more 
plausible here. 

By fixing the value of T^st to 13.8 K, i.e., to the average 
Jkin derived by Harju et al. (1993), we also computed the map 
of dust emissivity spectral index, /3 (Fig. [5}. The f3 varies from 
~ 1.3 to 2.8 inside the 3<x SABOCA contour. For comparison, 
through a multiwavelength study of the starless core TMC-1C, 
Schnee et al. (2010) derived the values 1.7 < /3 < 2.7, the 
most likely one being ~ 2.2. Shirley et al. (2011), by study- 
ing the Class source B335, found a relatively steep submm 
opacity power-law index of/? = (2.18 - 2.58)^3q. From Fig. [5] 
we see that f3 appears to rise to values > 2 towards those parts 
where T^st is increased. Of course, our present data does not 
allow to draw any firm conclusions and one must be careful 



not to overinterpret the current results, but it can still be spec- 
ulated that the increase of f3 is caused by the destruction of 
dust grains. This would be in line with the destructive outflow 
feedback discussed above. Depending on the velocity of the 
jet-induced shock, dust grains can be either partially or com- 
pletely destroyed by sputtering, grain-grain collisions, and shat- 
tering (e.g., Uones et al. 19941 IVan Loo et al. 20l"3l . Moreover, 
the processing by both UV photons and cosmic rays can cause 
non-thermal desorption of volatile species from the icy grain 



mantles into the gas phase (e.g., Hartquist & Williams 1990 



Roberts et al. 2007). However, it is also possible that higher val- 
ues of /? are actually caused by the growth of icy mantles on 
dust grains (e.g., [Schnee et al. 2010al and references therein), 
but this seems to contradict the chemical features discussed 
earlier. Conversely, observational results of small values of jS 
are often interpreted to be indicative of grain growth (e.g., 
ITesti et al. 20031 IDraine 20061 IKwon et al. 20091 . For example, 
Kwon et al. (2009) found that for their sample of Class sources 
(3 < 1, implying larger grain sizes due to gas accretion and/or 
coagulation. We note, however, that a visual inspection of the 
Spitzer/IRAC 3.6-pm image of 105399 reveals the presence of a 
coreshine-like emision (Fig.[3j. The 3.6-yum coreshine originates 
in dust-grain scattering of the background radiation, and is an in- 
dication of the presence of large, micron-size grains whose abun- 
dance is roughly given by the MRN (Mathis et al. 1977) grain- 
size distribution (Steinackeret al. 2010). For example, Stutz et 
al. (2010) found this towards the Bok globule CB 244 (coin- 
cidentally, it resembles our source in the sense that it contains a 
Class protostar and a starless core). Finally, the anti-correlation 
between Fdust and/? claimed in some observational studies might 
be artificial and caused by noise effects and Jdust-/? degeneracy 
instead of being a real physical trend (e.g., |Kelly et al. 2012| and 
references therein). 



5.2. IRAS 05399-01 21 /SMM 1 - a double core system 

One of the original aims of the present high-resolution SABOCA 
study was to search for substructures within 105399 and SMM 
1. The fact that SMM 1 remains a single source is consistent 
with the finding by Schnee et al. (2010b) that none of their 11 
starless cores break into smaller components at 5" resolution. 
Although, higher interferometric resolution may be necessary to 
see substructure at such scales (Off ner et al. 2012at . 

The I05399/SMM 1 -system appears to be a filamentary-type 
object not far from equilibrium (slightly thermally supercritical), 
and the origin of the two subcores, projectively separated by 0.14 
pc, seems to be caused by cylindrical Jeans-type fragmentation. 
We have not applied a correction factor for inclination, but the 
long length of the outflow driven by 105399, and its possible 
interaction with SMM 1 suggest that the system is aligned close 
to the plane of the sky. 

Although the fragmentation of the studied system can 
be explained in the context of a self-gravitating isothermal 
equilibrium filament, the estimated density profile n(r) oc 
r~( 2 - 3 ±o<>), albeit quite uncertain, is shallower than what would 
be theoretically expected, i.e., ~ r~ 4 (Ostriker 1964; see also 
lArzoumanian et al. 201 UlPineda et al. 201 11 1. Instead, it is close 
to the standard singular isothermal sphere profile of n(r) oc r~ 2 
(IShu 1 9771 1 . However, Fiege & Pudritz (2000) found that most 
of their models of isothermal filaments that are threaded by heli- 
cal magnetic fields and are in gravitational virial balance, can 
yield density profiles that decline like ~ r~ L8 to ~ r~ 2 (see 
also lLadaetal. 1999k Similarly, Caselli et al. (2002) found that 
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starless cores can be modelled with radial density profiles of 
n(r) oc r outside ~ 0.03 pc, whereas protostellar cores are 
better described with steeper single power-law density profiles 
(p > 2). For comparison, Shirley et al. (2002) modelled seven 
Class objects using a single power-law density profiles and 
found that the best-fit power-law index is p - 1.8 ± 0.1. The 
same density profile {p — 1 .8) was found by Schnee & Goodman 
(2005) for TMC-1C through SCUBA 450/850 pm observations, 
resembling our 350/870 pm analysis. 

If both 105399 and SMM 1 were formed simultaneously 
through fragmentation, it is interesting why they are at so dif- 
ferent stages of evolution: the former is already a relatively 
evolved protostar, while the latter one is still starless. This is 
not a unique case in this sense; a similar situation holds for 
the double-core globules CB 26 and 244 recently studied by 
Launhardt et al. (2013). Also, the filamentary clump associ- 
ated with the Class protostar IRAS 05405-0117 in Orion 
B9 shows starless substructure in the N2H + (1 - 0) map by 
Caselli & Myers (1994) and in submm dust emission maps 
(Papers I and III). Miettinen (2012) suggested that the core 
formation in Orion B9 could have been triggered by the feed- 
back from NGC 2024. Externally induced compression could 
have led to a phase of rapid accretion (Henne belle et al. 2003 . 



Motoyama & Yoshida 2003 1, which might be able to explain 



the advanced evolutionary stage of 105399, but it is unclear 
why SMM 1 would have not experienced the same fate. The 
length of the outflow driven by 105399 supports the hypothe- 
sis that the central source is a strong accretor. Because 105399 
lies at the tip of both the major parent filament (see Fig. 9 
in IMiettinen 2012l l and the smaller scale elongated structure 
studied here, it could represent a point where gravitational fo- 
cusing has enabled the rapid evolution into protostellar phase 
dTobin et al. 20101 >. Again, why this would have not happened 
for SMM 1 remains unclear. 

The virial mass we derived for SMM 1 in Paper II is M wn - - 
11.1 ± 0.3 M0, where it was assumed that n(r) oc r , If this 
is compared with the revised 870-^/m core mass derived here, 
we obtain the virial parameter a v { r = M^/M - 1.7 ± 0.3 
dBertoldi & McKee 1992b . If the density profile of SMM 1 is 
of the form n(r) oc r~ 23 , the above M vn - should be divided 
by a factor of 2.62, and, as a cigar-shaped prolate spheroid 
with the axial ratio of about 1.8, a further division by a fac- 
tor of ~ 0.8 should also be applied dBertoldi & McKee 19921 
ISchnee & Goo dman 2005). These correction factors yield the 
value Qf v ir = 0.8 ±0.2. SMM 1 is therefore very likely to be gravi- 
tationally bound (a v i r < 2) and near virial equilibrium (a v ir = 1). 
Also the high density and degree of deuterium fractionation in 
SMM 1 support its prestellar nature. The fact that 105399 has 
collapsed faster than SMM 1 could mean that the latter has 
stronger internal forces resisting gravity. One might also con- 
clude that the outflow from 105399 has triggered the formation 
of SMM 1, but that it has not collapsed to form a protostar yet. 



The source is filamentary in shape, and there is a protostel- 
lar core at the nortwestem tip (105399) and a starless core at the 
southeastern end (SMM 1). By comparing the total line mass of 
the SABOCA filament to its critical equilibrium value, we con- 
clude that as a whole it is slightly thermally supercritical. The 
observational fact that there appears to be two subcores projec- 
tively separated by 0.14 pc can be understood in terms of Jeans- 
type cylindrical fragmentation. A simple analysis suggests that 
the radial density profile is shallower than for an isothermal un- 
magnetised filament in equilibrium (~ r v 23 -m) vs. ~ r~ 4 ). 

The SED we constructed for 105399 suggests that it is near 
the transition phase between Stages and I, in agreement with 
our previous results. To our knowledge, we have presented the 
first Spitzer/IRAC images of the spectacular outflow driven by 
105399. The quadrupolar- type jet morphology seen in the 3.6- 
pm image supports the possibility that 105399 is actually a bi- 
nary protostar. 

If 105399 and SMM 1 were formed simultaneously through 
fragmentation of the parent filament, the distinct evolutionary 
stage between the two could be the result of rapid mass accre- 
tion onto 105399 in a high-pressure environment. In principle, 
this could be the result of core formation triggered by external 
feedback from NGC 2024 as suggested by Miettinen (2012), but 
why SMM 1 would not have experienced the same effect is un- 
clear. The temperature map we computed (under the assumption 
of constant dust opacity spectral index /3 - 1.8) shows a few 
warmer parts, one associated with 105399 and the other at the 
southeastern end. These could, in part, be imprints of the out- 
flow feedback from 105399, as supported by our earlier results on 
the chemical properties of the source (e.g., low CO depletion). 
However, it cannot be excluded that NGC 2024, some ~ 30' (~ 4 
pc) southwest of I05399/SMM 1, is proving external heating, al- 
though the source's temperature gradient does not support this 
idea. 

One of our original aims here was to search for substructures 
or dense condensations within the cores. Given the results pre- 
sented in Paper III (many of the Orion B9 cores are split into 
substructures consistent with Jeans instability), the high spatial 
resolution of the present SABOCA data, ~ 3 400 AU, could have 
revealed some further sub-fragmentation. However, the cores 
105399 and SMM 1 do not break up into sub-condensations but 
the system remains a double source at the current resolution. 
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6. Summary and conclusions 

We have mapped the dense core system IRAS 05399/SMM 1 
in Orion B9 at 350 and 870 pm using the bolometer cameras 
SABOCA and LABOCA on APEX. The present study demon- 
strates the complementarity of these two instruments, similarly 
to what has been done in other studies using SCUBA at 450 and 
850 pm on JCMT In particular, the present data allowed us to 
construct the maps of dust colour temperature and column den- 
sity at relatively high angular resolution (19'.'9). 
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Miettinen & Offner: SABOCA and LABOCA imaging of IRAS 05399-0121 
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Fig. 2. Top panels: Spitzer/MIPS 24- and 70-jum images of the I05399/SMM 1 -system. The images are shown with linear scaling, 
and the colour bars indicate the surface-brightness scale in MJy sr~' . The images are overlaid with contours of SABOCA 350-/im 
dust continuum emission as in Fig. Q] In the left panel, the yellow arrow starts from the base of the HH 92 jet driven by 105399 
(Bally et al. 2002 1. Bottom panel: Spitzer/IRAC three-colour composite image showing the I05399/SMM 1 -system in the centre 
and part of the outflow features associated with the HH 92 jet. The 3.6, 4.5, and 8.0 pm emission is coded in blue, green, and red, 
respectively, and the colours are shown in logarithmic scales. The contours show the 350-/im emission as in the top panels. A scale 
bar indicating the 0.5 pc projected length is indicated. Note that the total length of the outflow, 34', corresponds to 4.45 pc at 450 
pc. 
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Miettinen & Offner: SABOCA and LABOCA imaging of IRAS 05399-0121 
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Fig. 4. Distributions of dust colour temperature {left) and H2 column density {right). The dust emissivity index, fj, was fixed at 1.8. 
The colour-scale bars on the right indicate the units in K and cm 2 , respectively. The plus signs show the target positions of our 
previous molecular-line observations. The diamond symbol indicates the 24-/mi emission peak of 105399. The direction towards 
NGC 2024 is indicated by an arrow in the left panel. The annuli used to compute the density profile shown in Fig. [6] are shown 
by the green circles, while the black solid line in the right panel shows the location of the slice used to extract the central column 
density discussed in Sect. 4.5. The resolution of the images is 19V86 or about 0.043 pc. 
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